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Exclusive Leptoproduction of Real Photons on 
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Abstract: Polarisation asymmetries are measured for the hard exclusive leptoproduction 
of real photons from a longitudinally polarised hydrogen target. These asymmetries arise 
from the deeply virtual Compton scattering and Bethe-Heitler processes. From the data are 
extracted two asymmetries in the azimuthal distribution of produced real photons about 
the direction of the exchanged virtual photon: ^UL with respect to the target polarisation 
and AhL with respect to the product of the beam and target polarisations. Results for 
both asymmetries are compared to the predictions from a generalised parton distribution 
model. The s'mtp and cos(Oi;^) amplitudes observed respectively for the .Aul and All 
asymmetries are compatible with the sizeable predictions from the model. Unexpectedly, 
a sm{2(j)) modulation in the asymmetry with a magnitude similar to that of the sin 
modulation is observed. 
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1. Introduction 



Generalised Parton Distributions (GPDs) [Q, ^, ^ encompass the familiar parton distribu- 
tion functions and form factors within a unified description of the nucleon. The study of the 
concept of GPDs has shown that they may provide a way to investigate the contribution 
of quark orbital angular momentum to the spin of the nucleon . In a frame in which the 
proton moves quickly in the "longitudinal" direction, GPDs contain correlated information 
on the monodimensional distribution of quark momentum fractions in that direction with 
their two-dimensional spatial distributions in the transverse plane Q. Access to GPDs can 
be achieved through the measurement of cross-sections and asymmetries in the exclusive 
production of photons and mesons leaving an intact nucleon Q . There are four chiral-even 
quark GPDs for the nucleon at leading twist: H, E, H, and E. The H and E distributions 
are quark helicity-averaged distributions, whereas the H and E distributions involve quark 
helicity differences. The H and H distributions conserve nucleon helicity, while the E and 
E distributions are associated with a change in nucleon helicity. 

Appearing in the descriptions of exclusive leptoproduction of photons or mesons, GPDs 
are dependent on four kinematic variables: x, ^, t, and Q^. The variables x and ^ are 
the average and half the difference of the longitudinal parton momenta as fractions of 
the "infinite" nucleon momentum in the initial and final states respectively, in a frame in 
which the initial proton moves quickly. The average fraction x is not directly experimentally 

2 

accessible and the skewness variable ^ is related to the Bjorken scaling variable xb = as 
^ — ^2^^^) ill the Bjorken limit where the virtuality of the exchanged photon = —q^^ 
oo, while xb and t are fixed. Here, p is the four-momentum of the target nucleon and q is 
the four-momentum of the virtual photon. There is currently no consensus as to how to 
define ^ precisely in terms of experimental observables at finite Q^] the results are reported 
as projections in x-q. The Mandelstam t variable is defined as the squared momentum 
transfer to the target nucleon, i.e. t = {p — p')'^ where p' is the four-momentum of the 
recoiling nucleon. The dependence of GPDs on is implicit because GPDs are subject 
to quantum chromodynamics evolution with Q^, which has been calculated perturbatively 
to leading order and next-to-leading order in the strong coupling constant [||, |2|, ^, ^, |^ . 

The Deeply Virtual Compton Scattering (DVCS) process {Ip — )■ / 7p) is the sim- 
plest of those currently measurable that provide access to GPD-related information. At 
leading-order, a quark in the nucleon absorbs the virtual photon 7* (which mediates the 
electromagnetic interaction between the incident lepton / and the target nucleon p) radi- 
ated by the incident lepton, and itself radiates a detectable real photon. This combination 
of initial (Ip) and final (I'jp) states appears in two experimentally indistinguishable pro- 
cesses, the DVCS and Bethe-Heitler (BH) processes. Unlike in the DVCS process, in the 
BH process the real photon is radiated from the initial or the scattered lepton. 

In this paper, asymmetries in the distribution of photons produced from a hydrogen 
target by a positron beam are measured, where both the beam and target are polarised 
parallel or antiparallel to the beam direction. The extracted asymmetry amplitudes are 
compared to those predicted by a GPD model Q. 

The four-fold differential cross-section for exclusive leptoproduction of real photons. 
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neglecting target polarisation components transverse to the direction of the virtual photon, 
is given as 

do- XBc'^Irp ^^^^ 



dxBdQ2d|t|d(/. 32(27r)4 Q4 ^/r:p72 ' 

where e is the elementary charge, the angle (f) is defined as the azimuthal angle between 
the plane containing the directions of the incident and scattered lepton trajectories and 
the plane containing the trajectories of the virtual and real photons and e = 2xb^, in 
which M is the mass of the nucleon. The square of the scattering amplitude |rp can be 
written as 



|Tp= ItbhP + |tdvCsP + 7"BHTdvCS + '^bh'^dvcs • (1-2) 

Although the squared- DVCS amplitude ItdvcsP is small relative to the squared-BH ampli- 
tude |tbhP in the kinematic range of the HERMES experiment ||l^, information on GPDs 
can be accessed via the Interference term denoted I, which arises from the interference of 
the scattering amplitudes of the two processes . To leading order in the electromagnetic 
coupling constant aem, the Fourier expansion of these quantities for a charged lepton beam 
and a target nucleon, where both the beam and the target are longitudinally polarised, 
reads: 



TBH 



i^BH 



Yl ^™ cos(n0) + P,PiY, cos( 



.n=0 



ra=0 



(1.3) 



I TDVCS I = Kbycs' 



^ c?YCS cos(n6) + P, s?ycs ^-^ , 



+Pz 



/ J ^ra,unp 
n=0 
1 



unp 



n=0 



n=l 



(1.4) 



n=0 



cos(n0) + ^ sJi.unp sin(n(/ 



n=l 



^Pz 



Pe. Y cj, LP cos(n(/)) + ^ s^p sin( 



no 



n=0 



n=l 



(1.5) 



Here Vi{(j)) and V2{(t>) are lepton propagators of the BH process, P^ and respectively 
represent the longitudinal polarisation and the charge of the lepton beam in units of the 
elementary charge and Pz represents the longitudinal polarisation of the target. The 
subscript unp (LP) denotes coefficients for an unpolarised (longitudinally polarised) tar- 
get. The terms i^BHj -^DVCS and Ki are kinematic factors: -ftTBH = 1/ \X%i 



1 



-f^DVCS = and Ki = 1/ {xb yt), where y is the fraction of the beam energy carried by 

the virtual photon in the target rest frame. The lepton propagators and BH coefficients 
can be calculated in QED with the latter also having a dependence on Fi and 



■'n,{unp|LP) 



Fo, the Dirac and Pauli form factors of the nucleon. The Fourier coefficients c^y-"^ ,j and 

^' n,(unp|LP) 

sJ^y^p|Lp) arising from the squared-DVCS term relate to a bilinear combination of GPDs, 
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whereas the Fourier coefficients (unp|LP) (unp|LP) from the Interference term relate 

to a hnear combination of GPDs. 

Two asymmetries in the azimuthal distribution of real photons for positron scattering 
from a longitudinally polarised proton target are presented in this paper: one single-spin 
asymmetry ^UL arising from the longitudinal polarisation of the target averaged over all 
beam polarisation states, and one double-spin asymmetry All arising from the longitudinal 
polarisation of both beam and target. These can be written as 



ALL{(t>, eg) 



> eg) + a' 



. eg)] - [(7^^ (0, eg) + o-^^ (0, eg)] 



[a^^ {(j), eg) + a^^ (</., eg)] + [a^^ eg) + a^^ (0, eg)] 



(1.6) 



i^DVCS Yl ^°LP^ 



no 



ra=l 



n=l 



1 



Kbh cj?,unp cos(n0) - egKi ^ cj^ cos(n(/>) 



n=0 



n=0 



ra=0 



[a^^{(f),eg) + a^^{^,eg)] - [a'^^{(t),eg) + a^^{(l),eg)] 
c/), eg)] + [o-^^ (0, eg) + a~*-^ {(/), eg)] 



> eg) + CT^ 
1 



+ i^DVCS Y "^".unp COs(n</>) 
=0 

;i.7) 
;i.8) 



Yl "^".LP cos(n</>) + i^DVCS Yl ^nlW cos{n(l)) - Y "^".LP cos(n0) 



n=0 



n=0 



n=0 



1 

7'i(9i)P2W 



i^BH ^ c™np cos(n(/)) - e^Ki ^ cos(n0) 



n=0 



n=0 



+ i^DVCS Y "^".unp COs(n(^) 



n=0 



where a denotes the cross-section from Eq. 1.1 



(1.9) 

(•^) represents the beam helicity state 
parallel (anti-parallel) to the beam momentum, and (=^) represents the target polarisa- 
tion state parallel (anti-parallel) to the beam momentum. In the case of positron scattering, 
eg = +1. 



2. Experiment and Event Selection 



Data were collected in 1996 and 1997 with the HERMES spectrometer |12] using a longitu- 
dinally polarised 27.6 GeV positron beam incident on a longitudinally polarised hydrogen 
gas target ||l^ in the HERA lepton storage ring at DESY. The integrated luminosity of the 



data sample analysed in this paper is approximately 50 pb^^ [14|. The average beam and 
target polarisations for the data are presented in Table ^. A brief description of the event 
selection is given in the following paragraphs. More details can be found in Refs. [15, 1^, 17|. 
Exactly one charged track identified as a positron and one photon were required within 
the acceptance of the spectrometer. The kinematic requirements imposed on each event 
are 1 GeV^ < < lOGeV^, 0.03 < xb < 0.35, W > 3 GeV, and < 22 GeV, where W is 
the invariant mass of the initial 7*p state and is the energy of the virtual photon in the 
target rest frame as determined from measurements of the energies of the positron before 
and after scattering. 
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Year 


Luminosity 
[pb-i] 


Beam Pol 
Pe < 


arisation 
Pe > 


Target Pol 
Pz <0 


arisation 

> 


1996 


12.6 ± 1.0 




0.514 ± 0.017 


-0.759 ± 0.042 


0.759 ± 0.042 


1997 


37.3 ± 3.2 


-0.531 ± 0.018 


0.497 ± 0.017 


-0.850 ± 0.032 


0.850 ± 0.032 



Table 1: The integrated luminosities in pb ^ of the data sets, with average beam and target 
polarisations. 



Photons were identified by selecting a single signal cluster in the electromagnetic 
calorimeter with no associated track in the rest of the spectrometer. The cluster was 
required to register an energy deposition of at least 5 GeV in the calorimeter in order to 
reduce background, and to register at least 1 MeV in the preshower detector in order to 
improve the resolution of the energy measurement. 

The polar angle between the virtual and real photons was required to be between 
5mrad and 45mrad. The lower limit on this requirement was imposed in order to en- 
sure that the azimuthal angle (p remains well-defined within the finite resolution of the 
spectrometer - it excludes very few events from the data sample. The upper limit on 
the requirement was determined by Monte Carlo (MC) studies [18|, which indicated that 
the data set above this value is dominated by background from the Bethe-Heitler process 
producing resonant states of the proton and from semi-inclusive meson production. 

As the recoiling proton was not detected, events were selected by requiring that the 
squared missing-mass M^ = {q + p- q'f of the ep — )• 67 X reaction corresponded to the 
squared proton mass within the spectrometer resolution, where q' is the four-momentum 
of the real photon. An "exclusive region" in the missing-mass distribution was determined 
from MC simulations described in Ref. |]l9[ . Corrections were applied to account for shifts 
in the distributions between data samples for different years. 

The value of t was determined from the kinematics of the scattered lepton and the 
real photon under the assumption of ep — ?• ep7 in order to avoid using the measurement of 
the energy of the real photon, which was subject to the largest experimental uncertainty. 
The requirement —t < 0.7 GeV^ was imposed in the exclusive region to further reduce 
background contamination |18|. 



3. Experimental Extraction of Asymmetry Amplitudes 



This paper presents Fourier amplitudes of the asymmetries of Eqs. 1.6 and 1.8, rather than 
the related Fourier coefficients of the cross-section, which appear in Eqs. 1.5-1.5. The link 
between them is clarified in the following section. 

Section ^ describes the selection of an event yield A^. Its expectation value can be 
written as 

{MiP,, P., cP,ee)) = CiP,) 7]icl))auxj{cl),e,)[l + P,AvL{cP) + P,PzALL{cP) + PeALu{cP)] , 

(3.1) 

where C is the integrated luminosity, rj the detection efficiency and cjuu denotes the cross- 
section for an unpolarised beam and an unpolarised target. The beam helicity asymmetry 
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^LU is not considered in this paper since the dataset presented here is a subset of data 
previously analysed with respect to the beam helicity pO[| . In analogy to the decomposition 
of the cross-section in Eqs. 1.3-1.5, the asymmetries defined in Eqs. |1.6| and |1.8| can be 
decomposed as 



jCos(0(/>) 
UL 



(3.2) 



n=l 



(3.3) 



n=0 



The azimuthal asymmetry amplitudes A in Eqs. 3.2 and 3.3 are extracted from the data 
using the maximum likelihood fitting method ||l^, 21], with each amplitude containing a 
combination of the Fourier coefficients from Eqs. |1.3| -p!T5| with the exception of A'^^^'^'^\ 
This term has no physical meaning. It is included only as a test of the normalisation of 
the function and is expected to be zero. 

Previously published measurements with transversely polarised [ p!g| ] and unpolarised 
targets ||2^, 22 1 were made with both electron and positron beams. The beam-charge 
dependence of the contribution from the Interference term to the cross-section then allowed 
the separation of the squared- DVCS and Interference terms via charge difference and charge 
average asymmetries. However, the longitudinally polarised hydrogen data set at HERMES 
was taken solely with a positron beam, so the separation of squared-DVCS and Interference 
terms is not possible. 

are distinguished mainly by the Fourier 



The asymmetry amplitudes in Eqs. 3.2 and 



coefficients in the numerators of Eqs. 1.7 and (see the first two columns in Table |2|). 
The extracted amplitudes may also be influenced by the (/)-dependent lepton propagators 



and/or the other (/)-dependent terms in the denominators of Eqs. 1.7 and |1.9| . 

Correlations between the asymmetry amplitudes determined by the fit were found to 
be small, and likelihood ratio tests for higher-order terms show a null result. Therefore, 



only the terms shown in Eqs. |3.2| and |3.3| were fitted to the data, and the result of that fit 
is presented in this paper. 

The correspondences between the individual asymmetry amplitudes from Eqs. |3.2| and 
3.3 and the Fourier coefficients in the decomposition of the differential cross-section (within 
a kinematic factor) from Eqs. 1.3-|1.5|, which are interpretable within the GPD framework, 
are clarified in Table |2[ The relation of these Fourier coefficients to GPDs is encompassed 
in C-functions pT| , p3| , p^ , 25 1. These functions depend upon combinations of Compton 
Form Factors (CFFs), which are convolutions of GPDs with hard scattering kernels, and 
consequently have real and imaginary parts. Furthermore, the CFF information contained 
in the C-functions enters with various degrees of suppression. Higher twist terms are, in 
general, suppressed by powers of ^ compared to leading twist (twist-2). There are two C- 
functions that appear at twist-2 in the observables reported in this paper, C^p and C^p'^^. 
While the latter is a bilinear combination of CFFs and their complex conjugates, the former 
is written 
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Asymmetry 


Contributory Fourier- 


Power of i 

y 


Dominant CFF 


Twist 


Amplitude 


Coefficients 


Suppression 


Dependence 


Level 


< sin (f) 


*1,LP 
„DVCS 

^1 T P 


1 
2 


Im 

ImCpvcs 


2 
3 


.sm(2(^) 


„DVCS 

t>0 T p 


2 
2 


ImC^p 
ImCDycs 


3 
2 


.sin(30) 
^TTT 


si 

3 LP 


1 


Im LP 


2 


.cos(09i) 


'-O.LP 
^DVCS 
•^O.LP 


1 
1 


ReC[p 
ReCDVCS 


2 
2 


< cos <j) 


•-l.LP 
^DVCS 


1 
3 


ReC[p 
ReCDVCS 


2 
3 


,cos(20) 


''2,LP 


2 


ReC[p 


3 



Table 2: The correspondences between the asymmetry amplitudes extracted from the data set 
and the Compton form factor dependent Fourier coefficients of the differential cross-section. The 
subscript T refers to C-functions that involve gluon transvcrsity and arc further suppressed by 

as 

7T 



where the dominant summand is FiT-L at HERMES kinematic conditions. The C^p-function 
is the only C-function that is dominated by CFF T-L, and therefore GPD H. Consequently, 
the asymmetry amplitudes shown in Table § that provide access to C^p offer the best 
possibility to constrain H. 

Examination of Table || reveals three asymmetry amplitudes (^^'l^'^'^''' ^ul''''' ^ll 
that have leading-twist contributions from the CLp-function, via the Fourier coefficients 
(cq LP' '^1 LP' ^1 Lp) ™ numerators of Eqs. L7 and L£. The A'^^^'^'^^ amplitude receives 
an additional twist-2 contribution from the C^jy^^-function, stemming from the Cq lp^ 
Fourier coefficient. The A^^^^'^^ and ^ll*^ amplitudes also receive a contribution from 
BH coefficients, Cq^p and c^lp respectively. The tangled mix of contributions to these 
amplitudes increases the difficulty of extracting information related to CFFs and therefore 
GPDs from it. The ^ul"^ ^^"^ ^ll''^ asymmetry amplitudes each receive contributions 
at the twist-2 level from the CLp-function, with a twist-3 contribution from the C^p^^- 
function. The dominance of the Interference term over the squared-DVCS term, combined 
with the additional ^-suppression arising from the twist-3 nature of the C^p''^^ contribution 

to these asymmetry amplitudes, makes ^UL*^ ^^'^ ^ll"^ simplest of the asymmetry 
amplitudes considered in this paper from which to extract GPD-related information. More 
specifically, it is expected that these measurements can be used to constrain the values of 
the CFF Ti. The ^UL*^ amplitude is sensitive to the imaginary part of the CFF T-L, whereas 
the amplitude is sensitive to the real part of the same CFF. 
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The higher order Fourier components of the asymmetries receive twist-3 or gluon- 
transversity contributions, or a combination thereof; the A^^'^'^^ and A'^^^'^^^ both have 
twist-3 contributions from the imaginary and real parts of Cj^p respectively. The A^^'^'^^ 
amplitude also has a contribution from the gluon-transversity-dependent C^^p^ i'^M^ while 
the Ay^^'^'^ amplitude depends on a gluon-transversity function from the Interference term, 
CJp LP . Both gluon-transversity functions are suppressed by ^ compared to the quark 
leading-twist case and are expected to be very small. To date these gluon contributions 
have not been included in any model predictions for any target state. 



4. Background Correction and Systematic Uncertainties 

The extracted asymmetry amplitudes are corrected in each kinematic bin for background 
contributions from semi-inclusive and exclusive neutral meson production. The method is 
described in Ref. |19| | and the corrected asymmetry amplitude Acorr is given by 

, ^ /semi^scmi ~ /cxcl^excl (AT\ 

^corr — TZn ■ l^-^J 

-L /semi Jexcl 

The fractional contributions to the data yield from semi- in elusive neutral mesons /semi and. 
exclusive pions /exci '^^^ estimated from MC simulations [^] . The average semi- inclusive 
contribution is 3.1%, with the greatest variation being across the projection in xb, from 
1.1 % to 8.0 %. The exclusive neutral pion contribution is less than 0.7% in every kinematic 
bin. This estimate of the exclusive fraction is supported by calculations from another 



model [26 1, and a data search for corresponding events at HERMES [27|. 

The background asymmetry amplitude from the semi-inclusive process ^semi is ex- 
tracted from data by reconstructing neutral pions from the two-photon decay, requiring a 
fractional energy E^o/z^ > 0.8 and an invariant mass 0.10 GeV < M^^ < 0.17 GeV. The 
asymmetry amplitude from the exclusive background A^xd cannot be extracted from data 
due to the small yield of exclusive neutral pions. Thus they are assumed to be zero with 
an uncertainty of corresponding to one standard deviation from a uniform distri- 

bution in the range [—1,1]. Half the effect of the total background correction is assigned 
as a systematic uncertainty. The statistical uncertainties of the background fractions and 



asymmetry amplitudes appearing in Eq. 4A are propagated through to the final statistical 
uncertainty of the amplitude. 

The resulting corrected asymmetry amplitudes presented in this paper contain contri- 
butions from the elastic BH/DVCS processes and processes in which the proton is excited 
to a resonant state. An MC simulation using a parametrisation of the form factor for 
the resonance region from Ref. [|2^ is used to estimate the fractional contribution of this 
resonance process and the individual cross-sections for single-meson decay channels are 
calculated with the MAID2000 program [|2^. The contribution to the exclusive sample 
averages 13%, with the greatest variation of between 5.6% and 33.6% being across the t- 
range covered in this analysis. We assign no systematic uncertainty due to this contribution 
and no correction is applied; this contribution is considered to be part of the signal. 
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Amplitude 


A =b (5stat. i (5syst. 


Missing 
Mass 
Shift 


Background 


Detector/ 
Binning 
Effects 


< sin 


-0.073 ±0.032 ±0.007 


0.002 


0.006 


0.002 


,sin(2(/>) 


-0.107 ±0.032 ±0.008 


0.000 


0.002 


0.007 


.sin(3(/>) 


0.015 ±0.032 ±0.009 


0.000 


0.001 


0.009 


.COS(O0) 


0.122 ±0.044 ±0.004 


0.001 


0.003 


0.003 


yi COS 


-0.047 ±0.062 ±0.029 


0.002 


0.003 


0.028 


.cos(20) 


0.067 ±0.062 ±0.007 


0.004 


0.003 


0.004 



Table 3: Results for the asymmetry amplitudes extracted over the kinematic range covered by 
HERMES, with their statistical and systematic uncertainties. For the latter, the various systematic 
uncertainty contributions are also presented. Not included are scale uncertainties of 4.2% (5.3%) 
arising from the target (beam and target) polarisation measurements. 

The extracted asymmetries are also subject to systematic uncertainties arising from 
the combined effects of detector misalignment, acceptance, smearing and finite bin width, 
in addition to the background correction described previously. The systematic uncertainty 
originating from the combined contribution is estimated from a simulation of the spec- 
trometer using a MC generator based on GPD parametrisation detailed in Ref. |^]. The 
resultant uncertainty is denoted "Detector /Binning Effects" and is shown in Table |^. 

Finally, we assign a systematic uncertainty due to the year-dependent shift in the 
exclusive region of the missing-mass distribution. We choose to assign one quarter of the 
effect that the shift has on the extracted asymmetry amplitudes. 

No systematic uncertainty is assigned due to luminosity differences because the lumi- 
nosity does not depend on the target polarisation and beam polarisation dependent weights 
are assigned to each event in the extraction. Possible uncertainties arising from extra QED 
vertices are neglected as the effects have been estimated to be negligible for HERMES |^0| . 

The magnitudes of all the contributions discussed in this section are given in Table |3|. 

5. Results and Model Comparison 

Figures || and § respectively show the asymmetry amplitudes ^4^^'^"''^^ and inte- 
grated over the HERMES acceptance as well as projected across the kinematic variables 
— t, xb and Q^. All values are summarised in Table ^. It should be noted that, within 
the HERMES acceptance, (xb) and (Q^) are highly correlated (see Table ^) and it is not 
possible to disentangle the dependences of the asymmetries on these quantities. The mea- 
surements are also subject to scale uncertainties from the measurements of the beam and/or 
target polarisations. These scale uncertainties are given in the captions to the figures. The 
fractional contribution to the data set from resonant state production (i.e. "Reso. frac") 
is estimated from MC simulations and shown in the bottom row of Figs. |l] and |. It is not 
known how this contribution may affect the values of the extracted asymmetries and since 
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here it is not experimentally separable from the non-resonant data, it is treated as a part 
of the signal. 
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Figure 1: Amplitudes of the target-spin asymmetry ^ul sensitive to a combination of the In- 
terference and squared-DVCS terms, for positrons incident on longitudinally polarised protons, as 
projections in — xb, and . The leftmost column shows the asymmetry values when the extrac- 
tion is performed in a single bin across the entire kinematic range of the data set. The error bars 
(open red bands) show the statistical (systematic) uncertainties and the solid blue bands represent 
the predictions from the "VGG Regge" GPD model described in Refs. There is an addi- 

tional 4.2% scale uncertainty due to the precision of the measurement of the target polarisation. 
The fractional contributions from resonance production estimated from an MC model are presented 
in the bottom panel. 



All amplitudes presented correspond to Fourier coefficients described in Ref. [11] relat- 
ing to the twist-2 and twist-3 CFFs shown in Table with the caveat that this relationship 



may be complicated by various Cji terms in the denominators of Eqs. [1.7] and [iT 

The first harmonic of the .AuL; when the extraction is performed in a single bin from all 
kinematics, exhibits the value AyL*^ = — 0.073 it 0.032 (stat.) ±0.007 (syst.). The kinematic 
projections provide no evidence of strong dependences on — i, xb, or . This asymmetry 
amplitude receives a mixture of twist-2 and twist-3 contributions, as shown in Table |2|. 
The primary contributor is C^p, which is twist-2 and is expected to dominate the twist-3 
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Figure 2: Amplitudes of the double-spin asymmetry All sensitive to the Interference, squared- 
DVCS and squared-BH terms in the scattering amplitude, for polarised positrons incident on lon- 
gitudinally polarised protons, as projections in —t, xb, and Q^. The leftmost column shows the 
asymmetry values when the extraction is performed in a single bin across the entire kinematic range 
of the data set. The error bars (open red bands) show the statistical (systematic) uncertainties and 
the solid blue bands represent the theoretical predictions from the "VGG Regge" GPD model de- 
scribed in Ref. [||, 31 . There is an additional 5.3% scale uncertainty due to the precision of the 
measurement of the beam and target polarisations. The fractional contributions from resonance 
production estimated from an MC model are presented in the bottom panel. 



contribution from C^p'~'^. 

The ^uL^^"^^ amplitude has the unexpectedly large value A^^'^'^^ = —0.106 it 0.032 it 
0.008 when extracted from the integrated kinematic range of the data set. The projections 
across —t, xb and in Fig. || show no obvious features. This asymmetry amplitude is 
expected to receive a mixture of quark twist-3 and gluon twist-2 contributions, and as such 
could have been expected to be small in the HERMES kinematic range. 

The CLAS collaboration also published extractions |32| of ^ul"^ ^^"^ A^^'^'^^ although 
without projections in —t, and across the kinematic region covered by CLAS. In 
Table |, the signs of the CLAS results are suitable for comparison with HERMES results. 
The sign of the ^UL*^ (^ul''^''^'*) amplitude has been inverted once (twice) to account for 
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the different beam charge used and to bring the angular definitions used by CLAS into 
concordance with the Trento convention used by HERMES Q. The single bin extraction 
value evaluated by CLAS has approximately three times larger magnitude than that eval- 
uated by HERMES, whereas the A^^'^'^^ asymmetry amplitude at CLAS is consistent with 
zero in contrast to the non-zero value extracted by HERMES. It should be noted that 
CLAS measurements are taken with a different beam charge and at larger average values 
of XB and t than HERMES, resulting in a different sensitivity to CFFs Ti and Ti as can be 
deduced from Eq. 3^. 

The A^^^'^^ and consistency test A'^^^^'^^ amplitudes were found to be consistent with 
zero over the kinematic range of the HERMES experiment. The former receives contribu- 
tions from the leading-twist LP'^^^^^tion. 

The results from the extraction of the double-spin asymmetry are presented in Fig. |2|. 
The twist-2 amplitude A'^-^^^'^^ is found to be 0.115 it 0.044 it 0.004 when the extraction is 
performed in a single bin across the entire kinematic range of the data set. No dependences 
are observed in the data across projections in —t,XB, or Q^. This asymmetry amplitude 
receives contributions from the squared-DVCS and Interference terms in the scattering 
amplitude, as shown in Table |2[ However, it also receives a dominating contribution from 
the Cq LP Fourier coefficient as shown in Eq. L9. It is therefore expected to be positive and 
non-zero, as confirmed by the data. 

The first harmonic of the double-spin asymmetry ^ll''^ ~ —0.054 it 0.062 it 0.029 
is compatible with zero when the extraction is performed in a single bin across the entire 
kinematic range of the data set. It exhibits no dependences on —t, xb, or Q^. It is expected 
to receive contributions from twist-2 and twist-3 terms, as shown in Table ^, and therefore 
may be non-zero at HERMES kinematics. Its value is dominated by the cJ^lp term arising 
from the squared-BH term in the expansion of the scattering amplitude. 

The A^^^'^'^^ amplitude is also compatible with zero. The kinematic projections show 
no dependence on —t,XB, or Q^. This asymmetry amplitude receives only twist-3 contri- 
butions, and so is expected to be small at HERMES kinematics. Unlike the single-spin 
asymmetry, there are no previous experimental measurements with which to compare the 
extracted double-spin asymmetry. 

All asymmetry amplitudes in Figs. |l| and ^ are presented in comparison with calcula- 
tions 1 31 1, labelled "VGG Regge", from the GPD model described in Ref. [^. Predictions 
from this model have been compared previously with HERMES asymmetries with respect 
to beam helicity and beam charge ||2^ and for a transversely polarised target [jl^ with 
limited success for certain choices of parameters. It remains the only predictive model 
available for comparison with data taken on a polarised target. The model is an implemen- 
tation of the double distribution concept [|l|, ^ where the kernel of the double distribution 
contains a profile function that determines the dependence on ^, controlled by a parameter 
b [^]. The b profile parameter can be set in the range b G [l,oo), where the GPD is inde- 
pendent of ^ in the limit 6—7-00. In this model, the ^-dependence and the x-dependence 
are entangled, but the t-dependence is factorised. Here a Regge-inspired hypothesis for the 
t-dependence is used. 

The model is used to calculate the differential cross-sections for the electroproduction 
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of real photons for each beam and target polarisation state. These differential cross-sections 
are used to construct asymmetries with which to compare the data. The cross-sections and 
therefore the asymmetries contain no provision for the production of resonant states of the 
target proton. The width of the theory bands in Figs. |l| and ^ originates from the variation 
of the unknown b profile parameters. These free parameters are independently controllable 
for valence and sea quarks, and can be used as a fit parameter for the extraction of GPDs 
from hard leptoproduction data. The other parameters of the model are chosen from those 
that best describe data that were previously published by HERMES [^] . 

Figure [l| shows that the reported amplitude ^UL*^ -^UL single-spin asymmetry 

is well described by the model, which predicts the kinematic trend and the magnitude of 
this amplitude with reasonable accuracy. However, the relatively large amplitude A^^'^'^^ 
is not described by the model, which predicts that the amplitude should be small and of 
opposite sign. Although Ref. |32] alludes to a large exclusive pion background in the CLAS 
kinematic region, which could cause a large sin(2i?^) amplitude, as previously stated a data 
search at HERMES |27] showed no such contamination. The model ^ itself suggests that 
any background due to exclusive pions should be small and this is supported by a different 
set of calculations from a different model |26|. This disagreement between the data and the 
model for this amplitude is surprising and could be interpreted as an unexpectedly large 
contribution from the gluon twist-2 amplitudes, which are not included in the GPD model 
shown. 

Figure |2| shows that the predictions made by the model regarding the magnitude and 
trends of the amplitudes of the ^ll asymmetry mostly agree with the data. It describes 
the positive, slightly increasing trend for the A'^-j^^^'^^ amplitude observed in the data across 
all three variables. The prediction of the model for values of the ^ll"^ compatible with 
the data within the uncertainties of the extraction. The model predicts a small value for 
^LL^^*^^ across all kinematic projections which is compatible with zero within experimental 
uncertainty. 

The model is successful in predicting four of the five experimental asymmetries pre- 
sented in this paper. However, the same model failed to predict previously published 
HERMES results for the beam helicity asymmetry. A further caveat is that the model 
does not include any effect from the resonance contribution to the data set. 



6. Summary 

Data on the hard exclusive electroproduction of real photons from the 1996 and 1997 
years of operation of the HERMES experiment are analysed. Two asymmetries in the 
azimuthal distribution of leptoproduced photons from a longitudinally polarised proton 
target are presented. The sini;^ modulation ^UL*^ single-spin asymmetry ^ulC^^) is 

found to be small and negative, with a weak kinematic dependence. The asymmetries are 
compared to calculations from the only available GPD model ^] for these processes 
on a longitudinally polarised target. The asymmetry amplitude is broadly compatible 
with the theoretical predictions from the model. The A^'^'^'^ amplitude is found to be of 
opposite sign and larger than expected. The double-spin asymmetry ALhifp) is extracted 
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Kinematic 
Bin 


[GeV^] 


(IB> 


[GeV^] 


.sin,* 
± t^stat. i ^syst. 


.sin(2,#.) 
± -^stat. ± f^ayst. 


.si„(3,f) 
=1= t^stat. i f^syst. 


integrated 


0.115 


0.096 


2.459 


-0.073 ± 0.032 ± 0.007 


-0.106 ± 0.032 ± 0.008 


0.015 ± 0.032 ± 0.009 


0.00 < -£ < 0.06 
U.UD < —t <^ U.i4 
0.14 < -t < 0.30 

0.30 < -t < 0.70 


0.031 
0.094 
0.201 
0.408 


0.079 
0. 103 
0.110 
0.123 


1.982 
2.531 
2.883 
3.587 


-0.008 ± 0.051 ± 0.012 

A n o c -L- n A r" '7 -L- a ait 
— O.UOO ± U.00( ± O.Ul/ 

-0.169 ± 0.070 ± 0.007 
-0.138 ± 0.109 ± 0.017 


-0.060 ± 0.050 ± 0.007 

A 1 1 A -L- A A C A _1_ A A 1 O 

-0.154 ± 0.069 ± 0.008 
-0.191 ± 0.116 ± 0.021 


0.024 ± 0.049 ± 0.005 

A f^l a -L- A A C A _1_ A AAA 

— 0.U16 ± U.ooy ih u.uuy 
0.020 ± 0.069 ± 0.016 
0.014 ± 0.115 ± 0.017 


0.03 < IB < 0.07 
0.07 < IB < 0.10 
0.10 < IB < 0.15 
0.15 < IB < 0.35 


0.096 
0.099 
0.123 
0.188 


0.054 
0.084 
0.121 
0.198 


1.437 
2.115 
3.108 
4.934 


-0.003 ± 0.053 ± 0.008 
-0.134 ± 0.064 ± 0.008 
-0.039 ± 0.070 ± 0.007 
-0.195 ± 0.093 ± 0.018 


-0.045 ± 0.053 ± 0.004 
-0.228 ± 0.064 ± 0.007 
-0.051 ± 0.069 ± 0.015 
-0.056 ± 0.089 ± 0.052 


0.031 ± 0.053 ± 0.004 
-0.077 ± 0.060 ± 0.009 
0.080 ± 0.069 ± 0.013 
0.101 ± 0.090 ± 0.012 


1.0 < < 1.5 
1.5 < < 2.3 
2.3 < < 3.5 
3.5 < < 10 


0.085 
0.098 
0.123 
0.178 


0.056 
0.079 
0.108 
0.170 


1.236 
1.862 
2.829 
4.865 


-0.043 ± 0.059 ± 0.004 
-0.079 ± 0.060 ± 0.007 
-0.111 ± 0.068 ± 0.007 
-0.054 ± 0.071 ± 0.008 


-0.093 ± 0.059 ± 0.007 
-0.149 ± 0.061 ± 0.007 
-0.042 ± 0.067 ± 0.009 
-0.130 ± 0.074 ± 0.019 


0.016 ± 0.058 ± 0.006 
-0.036 ± 0.060 ± 0.007 
0.059 ± 0.066 ± 0.009 
0.034 ± 0.074 ± 0.010 


Kinematic 
Bin 


i-t) 
[GeV^] 


<IB> 


<Q'> 
[GeV^] 


.cos(0*) 
± ^stat. ± ^syst. 


.COB* 

± ^stat. ± (^syst. 


.cos(2*) 
± ^stat. ± <^syst. 


integrated 


0.115 


0.096 


2.459 


0.115 ± 0.044 ± 0.004 


-0.054 ± 0.062 ± 0.029 


0.095 ± 0.062 ± 0.007 


0.00 < -t < 0.06 
0.06 < — i < 0.14 
0.14 < -t < 0.30 
0.30 < -t < 0.70 


0.031 
0.094 
0.201 
0.408 


0.079 
0.103 
0.110 
0.123 


1.982 
2.531 
2.883 
3.587 


0.129 ± 0.068 ± 0.010 
0.197 ± 0.080 ± 0.007 
-0.113 ± 0.095 ± 0.008 
0.237 ± 0.162 ± 0.009 


-0.012 ± 0.094 ± 0.010 
—0.021 ± 0.112 ± 0.022 
-0.179 ± 0.137 ± 0.044 
-0.065 ± 0.235 ± 0.079 


0.104 ± 0.097 ± 0.006 
0.031 it 0.114 ± 0.014 
0.206 ± 0.136 ± 0.007 
0.020 ± 0.211 ± 0.009 


0.03 < IB < 0.07 
0.07 < IB < 0.10 
0.10 < IB < 0.15 
0.15 < IB < 0.35 


0.096 
0.099 
0.123 
0.188 


0.054 
0.084 
0.121 
0.198 


1.437 
2.115 
3.108 
4.934 


0.137 ± 0.076 ± 0.014 
-0.111 ± 0.083 ± 0.012 
0.265 ± 0.095 ± 0.020 
0.192 ± 0.125 ± 0.070 


-0.108 ± 0.108 ± 0.013 
-0.023 ± 0.123 ± 0.018 
-0.169 ± 0.135 ± 0.016 
0.279 ± 0.178 ± 0.047 


0.094 ± 0.102 ± 0.003 
0.171 ± 0.117 ± 0.008 
0.087 ± 0.132 ± 0.012 
-0.193 ± 0.176 ± 0.017 


1.0 < < 1.5 
1.5 < < 2.3 
2.3 < < 3.5 
3.5 < < 10 


0.085 
0.098 
0.123 
0.178 


0.056 
0.079 
0.108 
0.170 


1.236 
1.862 
2.829 
4.865 


0.118 ± 0.080 ± 0.009 
0.061 ± 0.082 ± 0.005 
0.126 ± 0.091 ± 0.008 
0.164 ± 0.102 ± 0.010 


-0.120 ± 0.114 ± 0.029 
0.015 ± 0.120 ± 0.024 
-0.092 ± 0.130 ± 0.023 
-0.019 ± 0.149 ± 0.016 


0.111 ± 0.115 ± 0.004 
0.196 ± 0.113 ± 0.006 
0.040 ± 0.130 ± 0.005 
-0.076 ± 0.149 ± 0.013 



Table 4: Results of the ^ul"^' ^ul^^'^'' ^'^'^ A^^'^'^'^ amplitudes of the single-spin asymmetry (top) 
and ^ll"^' ^'^'^ ^^'^(2'^) amplitudes of the double-spin asymmetry (bottom) with statistical 

and systematic uncertainties and average kinematics for polarised hydrogen data. These are shown 
integrated over the HERMES acceptance and for each ~t, xb, and bin. The correlation between 
(xb) and (Q^) can be observed in the average kinematic values for each bin. There is an additional 
scale uncertainty of 4.2% (5.3%) affecting the ^™("'^) [A'^-^^"''^^ ) amplitudes due to the uncertainty 
in the measurement of the target (beam and target) polarisations. 



Experiment 


[GeV2] 




m 

[GeV^] 


^sin</> _|_ ^^^^^ _|_ Jgyg^^ 


<sin(2</.) , . , . 
^UL -■- "stat. ± fsyst. 


HERMES 


0.12 


0.10 


2.46 


-0.073 ± 0.032 ± 0.007 


-0.106 ±0.032 ±0.008 


CLAS 


0.31 


0.28 


1.82 


-0.252 ±0.042 ±0.020 


-0.022 ±0.045 ±0.021 



Table 5: The single-spin asymmetry defined in Eq. |l.6| , extracted at HERMES and CLAS with 
the average kinematic values for the respective datasets. The sign of the CLAS sm4> result has 
been inverted to become consistent with the different angular definitions used by HERMES and the 
different beam charge used in the two measurements, as described in the text. 



for the first time. It is found to be consistent with the small asymmetry values predicted 
by the GPD model. Both of the asymmetries have dominant contributions from the real 
or imaginary parts of the CFF T-L^ and thus provide potentially the best access to the GPD 
H at HERMES kinematics. 
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